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O . ABSTRACT 

Q^,) ' Using hybrid N- body /semi- analytic substructure simulations, we have devised a 

method to select galaxies that are isolated in their dark matter halo (N = 1 systems) 
and galaxies that reside in a group of exactly two (TV = 2 systems). Our N = 2 sys- 
tems are widely-separated with comoving projected separations up to ~200/i _1 kpc, 
where close galaxy-galaxy interactions are not the dominant physical process at work. 
We apply our selection criteria to two volume-limited sample of galaxies from the 
NYU Value-Added Galaxy Catalogue based on SDSS DR6 with M r -51og 1Q h < -19 
and —20. The effects of a simple change in environment on star formation properties 
are isolated by correcting the populations for contamination by systems in other en- 
vironments and by eliminating differences in stellar mass. For satellite galaxies in a 
!> ■ group of two, we find a red excess of 0.15 ± 0.01 and 0.14± 0.01 for the —19 and —20 

samples, respectively, relative to isolated galaxies of the same stellar mass. A simple 
galaxy formation model that assumes (1) N = 1 galaxies are the progenitors of N = 2 
galaxies in a statistical sense and (2) cold gas is instantaneously and completely re- 
moved from an infalling satellite galaxy immediately upon accretion into another halo 
predicts a red excess of 0.62 ±0.01 and 0.54 ±0.01 for the —19 and —20 samples, re- 
spectively. Thus, an immediate-rapid star formation quenching scenario is inconsistent 
with these observations. A delayed-then-rapid star formation quenching scenario with 
a delay time of 3.3 and 3.7 Gyr for the —19 and —20 samples, respectively, yields a 
red excess prediction that agrees with the observations. The observations also reveal 
that central N = 2 populations have a slight blue excess of 0.06 ± 0.02 and 0.02 ± 0.01 
rS • for the —19 and —20 samples, respectively, relative to N = 1 populations of the same 

^ | stellar mass. Additionally, we find red satellites galaxies slightly tend to pair with red 

centrals while blue satellites slightly tend to pair with blue centrals demonstrating 
galactic conformity is somewhat present even in groups of two. Taken together, our 
results demonstrate that environmental processes strongly influence even the sparest 
groups of luminous galaxies and in-depth studies of these simple systems are an es- 
sential step towards understanding galaxy evolution in general. 

Key words: galaxies: evolution - galaxies: haloes - galaxies: interactions - galaxies: 
clusters: general - galaxies: statistics -methods: statistical 



OO 

in 



1 INTRODUCTION 

The overall galaxy population out to z ~ 1 divides i nto two 
distinct types: early-type and late-type galaxies ( e .g. [H ubble 
19261 : IStrateva et alj I2OOI): felanton et all 12003; iBell et al. 
20041 : iTanaka et al.ll2005l : iBaldrv et al.ll2006l : ICooper et al. 

* E-mail: c.trinh@physics.usyd.edu.au 



2006). Although the exact definition of the two types varies 
in the literature, late-type galaxies tend to have spiral mor- 
phology and are blue due to a highly active recent star for- 
mation history. On the other hand, early-type galaxies tend 
to have elliptical morphology and are red due to a lack of 
recent star formation. 

From observations of galaxies at different redshifts we 
know that the stellar mass on the red sequence has increased 
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by a factor of 2 since z ~ 1 (|Bell et al.ll20ol ). Along with 
numerous other observations this suggests that physical pro- 
cesses have been at work since z ~ 1 transforming late-type 
galaxies into early-type galaxies. Young clusters at z ~ 0.5 
contain many spiral galaxies and very few SO galaxies while 
local clusters have a ratio of SO to elliptical galaxies greater 
by a factor of 5 and a similar dec rease in the number of spi- 
ral galaxies |Dressler et al.|[l997l ). Also, the galaxy luminos- 
ity function density normalisation for blue galaxies has re- 
mained essentially constant, while it has double d (~ 0.5 dex) 
for red galaxies since z ~ 1 jFaber et al.ll2007l ). 

In the current A Cold Dark Matter (ACD M) picture 
of structure formation (|Blumenthal et aTlll984r i. dark mat- 
ter collapses into haloes, trapping baryonic matter at early 
epochs. In galaxy-size dark matter haloes, cool baryonic ma- 
terial falls in while retaining angular momentum to form 
luminous spiral galaxies. Galaxy groups and cluster form 
hierarchically when smaller dark m atter haloes combine 
to form larger dark matter haloes JWhite fc Reed Il978l : 
iKauffmann et al.lll993l ; ICole et al]|200Ch . This picture mo- 
tivates a number of physical processes that may alter a 
galaxy's morphology and/or quench star formation through 
the removal and/or depletion of the galaxy's gas reservoirs 
bring about a late-type to early-type transition. 

The merger of two dark matter haloes will eventually 
result in the mer ger of their central gala xies as a result of dy- 
namical friction l|Chandrasekharl [l943) , potentially perma- 
nently altering both morphology and star formation rates. 
Even before the central galaxies merge they are altered by 
interactions, especially if a smaller dark matter halo has 
fallen into a large group or cluster-size halo with a dense intr- 
acluster medium and many other un-merged satellite galax- 
ies. These environmental interactions broadly fall into two 
categories: gravitational tides and gas collisions. 

Tidal interactions due to the main cluster potential on 
the infalling satellite galaxy may strip it of stars and/or 
gas leading to changes in morphology and/or star forma- 
tion rates |Fuiitdfl99sh . Tidal interactions from a close 
galaxy pass also triggers inflows of gas into the centres of 
galax ies, inducing short-lived starburst (|Mihos fc Hernquistl 
1 19961 ') . Once the galaxy has exhausted its supply of gas, 
it will quickly redden unless it can replenish its supply 
of gas. One mechanism that can prevent the accretion of 
new gas is feedback from an active galactic nucleus (AGN, 
e.g. iDi Matteo et alj|2005 iBower et ai1l2006l ; ICroton et ail 
2006). The energy released by an AGN generates an outflow 
of gas preventing accretion, which may transform a blue, 
active galaxy into a red and dead galaxy although the de- 
tails remain unclear. In addition to tidal interactions with 
the main cluster potential, a satellite galaxy will tidally in- 
teract with the other un-merged satellites within the halo. 
All these close high-speed encounters, which occur approx- 
imat ely once per Gyr, are ref e rred to as galaxy h arassment 
(e.g iFarouki fc Shapiro! Il98ll ; iMoore et all 1 19961 ) and can 
transform small disc galaxies into dwarf elliptical or dwarf 
spheroidal galaxies. 

The collisional interaction of the hot gaseous intraclus- 
ter medium of the larger host halo with the cold disc gas 
and/or the hot halo gas reservoir of the smaller halo's central 
galaxy may rapidly strip it of either g as reservoir in a pro - 
cess known as ram-pressure stripping (|Gunn fc Gotdfl972T ). 
High-resolution hydrodynamical simulations have shown 



that ram-pressure stripping can remove the entire H I gas 
content within 100 Myr from a luminous spiral galaxy like 
the Milky Way leading to red der colours from a lack of star 
formation (Q uilis et al.ll20o"ol ). However, the morphology of 
the original stellar disc is unaffected by ram-pressure strip- 
ping, which results in the spiral galaxy being transformed 
into an SO-like galaxy. If only the hot halo gas reservoir of 
a satellite galaxy, which replenishes the cold disc gas that is 
converted into stars, is removed by ram-pressure stripping 
or tides, the result is a gradual decline of star formation be- 
cause it continues until the cold disc gas is exhausted. To 
distinguish it from the removal of cold disc gas, which results 
in a rapid truncation in s tar formation, this is usually re- 
ferred to as strangulation (|Larson et al.lll980l ; iBalogh et al.l 
2000). Strangulation should result in a similar transforma- 
tion of a spiral into an SO-like galaxy, except on a different 
time-scale than ram-pressure stripping of cold disc gas. 

If an infalling satellite galaxy is able to withstand the 
environmental interactions discussed above, a major or mi- 
nor m erger with another satellite galaxy (M akino fc Hutl 
Il997h or the central galaxy of the host halo may result 
in a merger remnant with drastically different morphology 
and/or star formation rate than either of the two progeni- 
tors. The severity of the change is highly dependent on the 
mass ratio of the two progenitors. If a spiral galaxy is in- 
volved in a minor merger, a temporary starbursts may be 
triggered from the compression of colliding cold gas and 
the disc morphology may be affected. However, morphol- 
ogy is definitely altered when two spiral galaxies with sim- 
ilar masses experience a major merger. In this case, the 
mer ger remnant may have sp heroidal or elliptical morphol- 
ogy ( Toomre fc Toomrdl 19721 ). If the gas reservoirs depleted 
in the triggered starburst are not replenished because of say 
AGN feedback, the major merger will have transformed two 
active spiral galaxies into a red and dead spheroidal or el- 
liptical galaxy. 

A fundamental goal of astrophysics is a compre- 
hensive understanding of the role mergers and envi- 
ronmental processes play in building and evolving the 
diverse set of galaxies that exist in the Universe. 
Studies of galaxy evolution span a range of environ- 
ments including iso l ated galaxies (e.g. lAllam et all 120051 ; 
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2012 
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12004 Gerke et al. 2005; 


Weinmann et al. 2006a b 
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Kimm et all 20091; [Skibba 


2009 
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iTanaka et al.ll2005l ; Ivon der Linden et al.ll2007l , 120101 ) . 

However, a general understanding of galaxy evolution 
requires an understanding of the simplest galactic environ- 
ments, which are isolated galaxies residing alone in their 
dark matter halo (N — 1 system) and systems of two galax- 
ies sharing a dark matter halo (N = 2 system), which 
includes galaxy pairs and groups of two galaxies. Isolated 
galaxies are the least ambiguous and are the most controlled 
environments in which to study galaxy evolution. In general, 
more isolated galaxies are more likely to have later-type mor- 
phologies, higher star fo r mation rates and blue r colours (e.g. 
iPostman fc Gellerlll984l ; iBlanton et a"i]|2005bl) . 
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In this paper, we present an preliminary study of the 
difference between the properties of isolated galaxies and 
groups of two galaxies, the simplest group environment 
and the simplest change in galactic environment. We use 
a cosmological model consisting of a hybrid iV-body/semi- 
analytic substructure simulation to understand and remove 
the contamination from galaxies in other environments al- 
lowing us to study the full, uncontaminated distributions of 
star-forming and morphological parameters instead of just 
the average trend. This approach differs from pr evious stud- 
ies t hat utilise group-finding algorithms (e.g. lYang et al.l 
2005), which have difficulty accountin g for interloper sys - 
tems. This technique has been used bv lBarton et al.l (|2007T ) 
to study close galaxy pairs, where galaxy-galaxy inter- 
actions are dominant. The dominant environmental pro- 
cess (es) at work in groups is sti ll an open question and a 
very a ctive field of research (e . g. IWeinmann et aL 2006a R 



van den Bosch et al. 
Pasauali et al.ll2010l ; 



2009; Kan g fc van den Boschl 120081: 
200sl; iKimm et al.|l2009l; ISkibball2009l; 

McGee et al.ll2009l . Km! ; IWetzel etalj|2012f ). The low-mass 
haloes of groups of two galaxies examined in this paper rep- 
resent the minimum mass where t he dominant group envi- 
ronmental process begins to occur l|McGee et al J [2009 V In- 
vestigations of groups of two aim to shed light on the subject 
and they afford a number of other advantages including the 
fact that ACDM makes robust p redictions for the m erger 
histories of N = 2 systems (e.g. IStewart et al.l [20081 ) and 
such systems are routinely s tudied in high-re solution hydro- 
dynamical simulations (e.g. ICox et aT1l2006h allowing for a 
direct comparison yielding insights assumptions and treat- 
ment of interactions, star formation and feedback. As such, 
groups of two galaxies are a valuable vehicle for the study of 
galaxy evolution and will contribute toward a comprehensive 
understanding of the topic in general. 

The layout of this paper is as follows. In Section [2] we 
discuss the cosmological model. Section [3] describes the se- 
lection of populations of isolated galaxies and groups of two 
galaxies starting with an analysis of the environments of the 
simulated galaxies in our cosmological model. Section 13.11 
contains a description of the observational data set result- 
ing from the application of our selection criteria to galaxie s 
from the Sloan Digital Sky Survey fSDSS. lYork et alj|2000h . 
Our method for correcting for the contamination by galax- 
ies in other environments is explained in Section 13.21 We 
look at the stellar mass distribution of isolated galaxies and 
groups of two galaxies and discuss our Monte Carlo mass 
resampling technique for producing populations of the same 
stellar mass in Section [3.31 Results are discussed in Section 
[4] In Section |4T] we concern ourselves with the difference in 
colour between isolated galaxies and groups of two galaxies. 
In Section 14.21 we discuss the red excess of satellite galax- 
ies in a group of two and its implications on star formation 
quenching time-scales. The blue excess of central galaxies in 
a group of two and its origins are explored in Section 14.31 
followed by a conclusion in Section [5] 



2 COSMOLOGICAL MODEL 

The main components of our cosmological model include a 
large A-body simulation of cold dark matter halo formation 
and a semi-analytic cold dark matter substructure model 



jZentner fc Bullock! liobll ; IZentner et ail 120051). This m odel 
has been previously discussed in Berrier et al.l |2006l ) and 
iBarton et al.l (|2007l ) and we provide a brief review here. 

The model uses A-body simulations to characterise 
the spatial and mass distributions of "host" dark matter 
haloes. By definition, the centres of host haloes do not lie 
within the virial radius of other haloes. Our TV-body sim- 
ulation follows the evolution of 512 dark matter particles 
until z = in a comoving box of volume 1 20 ft" 1 Mpc 3 us- 
ing th e Adaptive Refinement Tree code of iKravtsov et all 
i|l997l ) in a standard ACDM cosmology with f2 ra = 0.3, 
Q.A = 0.7, ft = 0.7, and erg = 0.9. The implied particle mass 
is m p ~ 1.1 x 10 9 ft -1 Mq and the simulation grid was re- 
fined down to a minimum cell size of ft pGa k — 1.8 h~ kpc 
on a side. For more de t ails re g arding the A-body simu- 
lation see IZentner et all (|2005l ). lAllgood et all (|2006l ) and 
IWechsler et all (|2006l ). 

A-body simulations predict the substructure content 
of host haloes, but suffer from various numerical resolu- 
tion limitations. Thus, we populate each host dark mat- 
ter halo in the simulation volume wi th subhaloes u si ng th e 
semi-analytic substructure model of IZentner et al.l (2005), 
which effectively has infinite resolution. We use the semi- 
analytic model randomly to generate four independently re- 
ali sed mass accretion h is tories using the stochastic method 
of ISomerville fc KolattJ l|l999h for each host halo of mass 
M at redshift z in our simulation volume. Then, we de- 
termine the orbital evolution for each infalling subhalo in 
the potential of the host from accretion until 2 = 0, dur- 
ing which it loses mass and its maximum circular velocity 
decreases as its profile is heated by interactions. If the max- 
imum circular velocity of a subhalo falls below 60kms _1 at 
any time, it is removed to mimic the dissolution of the ob- 
servable galaxy as a result of these interactions. Each mass 
accretion history results in a different subhalo population 
in the simulation volume. We interpret these realisations as 
four independent cosmological volumes with identical large- 
scale structure, but different small-scale structure. All four 
realisations will be used in the statistical analysis that fol- 
lows. 

We assume each subhalo in the simulation volume con- 
tains a luminous satellite galaxy and each host halo has a 
luminous central galaxy. For host haloes, we use the max- 
imum circular velocity as a proxy for luminosity and as- 
sume a monotonic relationship between the two. For sub- 
haloes, w e use the maximu m circular velocity at accre- 
tion (sec Berr ier et alj [2006). In other words, all haloes 
above a given cutoff Knaxwill correspond to a popula- 
tion of galaxies brighter than some absolute magnitude. 
By exploring the number density of haloes in the model 
as a function of Vmax, we find that the number density 
of haloes with Vmax > 140 (196) kms -1 matches the ob- 
serve d number density of SDSS Da ta Release 6 (SDSS 
DR6, Adc lman- McCarthy et all 120081 ) volume-limited sam- 
ples with M r - 5 log 10 ft < -19 (-20). 

With this cosmological model we construct mock red- 
shift surveys of galaxies from SDSS DR6. With these mock 
volume-limited redshift surveys we quantify the number of 
luminous galaxies that occupy every host halo in the simu- 
lation, N, an unobservable parameter that is a key measure 
of the environment of a galaxy. We analyse the correlation 
between N and observable measures of environment, such as 
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Figure 1. (Left) Fraction of simulated galaxies in our mock volume-limited redshift survey (M r -51og 10 h < -20 or V mBX > 196 kms" 1 ) 
that reside in a dark matter halo of a given multiplicity TV, as a function of the comoving projected distance to their nearest neighbour, 
Dff. Each bin in _Djv is normalised to 1. Djv alone cannot be used to define samples of galaxies in N = 2 systems from volume-limited 
redshift surveys. Galaxies with small Djq values are predominantly in haloes with a total of 9 or more galaxies. (Right) Fraction of 
simulated galaxies with exactly one neighbour within 700 kpc (/V700 = 1) that reside in a dark matter halo of a given multiplicity N, 
as a function of -Djv- Very few galaxies with -/V700 = 1 reside in haloes with more than two galaxies. Galaxies with -/V700 = 1 and large 
values of Djv have a very high probability of being isolated in their halo and those with -/V700 = 1 and small values of Djv have a high 
probability of being in an isolated galaxy pair. 



projected nearest neighbour distance. This analysis reveals 
selection criteria for defining relatively pure populations of 
isolated galaxies (N — 1 systems) and groups of two galaxies 
(N — 2 systems) from volume-limited redshift surveys and 
gives statistical information on other unobservables of these 
populations. 



3 SELECTION CRITERIA FOR N = 1 AND 
N = 2 SYSTEMS 

First, we analyse the environment of galaxies in our mock 
volume-limited redshift surveys in order to define selection 
criteria for isolated galaxies (TV = 1) and groups of tw o 
galaxies (N = 2) in SDSS. Following Barton et all |2007l ). 
for each "galaxy" (halo or subhalo with V max > 140 or 
196 km s -1 ) in our mock volume-limited sample, we com- 
pute Djv, the comoving projected distance to the galaxy's 
nearest neighbour within AV < 1000 km s -1 , using peri- 
odic boundary conditions to fully sample its environment 
when necessary. We also compute N700 , the total number of 
galaxies within a comoving projected distance of 700 ft -1 kpc 
within AV < lOOOkms -1 . Lastly, we measure N, the total 
number of galaxies that lie within the same host dark mat- 
ter halo, which is unobservable in galaxy redshift surveys. 
Constraints on the observables Djv and N700 can identify 
galaxies in redshift surveys with a high probability of being 
in an N = 1 or N — 2 system, by applying them simultane- 
ously to the d ata and the simulatio ns. 

Following iBarton et ail (|2007m . but with a different 
Knax cutoff, Figure [T] (left) shows the fraction of galaxies 
with V max > 196 km s -1 in our mock volume-limited sample 
that reside in a halo of a given multiplicity iV as a function 
of Djv. Galaxies that are relatively far from their nearest 



neighbour are overwhelmingly in N = 1 systems. Galaxies 
that are relatively close to their nearest neighbour are pre- 
dominantly in haloes with nine or more galaxies in the same 
dark matter halo. The results are qualitatively similar for 
Knax> 140 km s" 1 . Thus, using Djv alone does not provide 
a sufficient means to select for N — 2 systems in volume- 
limited redshift surveys. 

However, these systems can be selected using Djv and 
-^700 together. FigurefJ] (right) shows the fraction of galaxies 
in our mock volume-limited sample with -/V700 = 1 that re- 
side in a halo of a given multiplicity TV as a function of Djv. 
The figure shows that very few galaxies with -/V700 = 1 re- 
side in haloes with N > 2. Moreover, galaxies with N700 = 1 
and large values of Djv are in N = 1 haloes and those 
with small values of Djv reside almost exclusively in N = 2 
haloes. Thus, we identify a population of galaxies with a 
high probability of residing in an N — 2 system can be de- 
fined by selecting galaxies from volume-limited redshift sur- 
veys with exactly one neighbour within a projected distance 
of 700 h^ 1 kpc and Djv smaller than some maximum value. 
The purity and contamination by interlopers in the N — 2 
sample may be determined from our cosmological model by 
quantifying the number of galaxies residing in N = 1, N = 2, 
and N > 2 haloes for a given maximum Djv. 

We choose to focus on groups of two galaxies because 
we are interested in studying the effects of groups processes 
on galaxy properties in the simplest group environment; by 
choosing a large maximum Djv, our N = 2 sample will con- 
tain few close pairs that are likely to be strongly interact- 
ing. In fact, the cosmological model shows that for Djv < 
200 h~ kpc, the fraction of iV = 2 galaxies that have had 
a close pass within 50 or 30 h~ l kpc within the last 0.5 Gyr 
is ~ 11 and 4 per cent for the —19 sample and 9 and 3 
per cent for the —20 sample. Table [1] lists the fraction of 
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Table 1. Sample definitions, sizes and purity. Symbol definitions are as follows: Dn - comoving projected distance to the nearest 
neighbour, TV700 - number of neighbours within a comoving projected distance of 700 /i — 1 kpc, -/Vgrjgg - number of galaxies in 
the SDSS sample, N s - lm - number of simulated galaxies in the cosmological model, Jn—1 - fraction of simulated galaxies in 
N = 1 systems, Jn=2 - fraction of simulated galaxies in N = 2 systems, /jv>2 - number of simulated galaxies in N > 2 systems 



Sample 


Magnitude Limit 


D N (ft- 1 kpc) 


N700 






/jv=i 


In=2 


/iV>2 


N = 1 


< -20 


> 400 


< 1 


26,192 


18,030 


0.995 


0.001 


0.004 


N = 1 


< -19 


> 400 


< 1 


7,266 


37,309 


0.999 


0.000 


0.000 


N = 2 


< -20 


< 250 


1 


3,229 


2,675 


0.203 


0.788 


0.008 


N = 2 


< -20 


< 200 


1 


2,549 


2,276 


0.144 


0.848 


0.007 


N = 2 


< -20 


< 150 


1 


1,830 


1,807 


0.106 


0.887 


0.007 


N = 2 


< -19 


< 250 


1 


1,096 


5,634 


0.369 


0.629 


0.002 


N = 2 


< -19 


< 200 


1 


918 


4,786 


0.289 


0.709 


0.002 


N = 2 


< -19 


< 150 


1 


721 


3,769 


0.199 


0.799 


0.002 



N = 1 haloes, /jv=i, the fraction of TV = 2 haloes, /jv=2, 
and the fraction of TV > 2 haloes, /jv>2, for a maximum 
D N of 250, 200, and 150 h' 1 kpc for both luminosities. The 
Dn < 200 hT 1 kpc population is our fiducial TV = 2 sample 
but we verify all results for the other two values. 

A population of galaxies with a very high probability 
of being isolated in TV = 1 systems results from selecting 
galaxies with TV700 < 1 and Dn > 400 h~ x kpc. The purity 
and contamination for this population for both luminosities 
are also listed in Table [T] Note that the N = 1 population 
defined in this way is very pure (99.5-99.9 per cent). 

3.1 TV = 1 and TV = 2 systems in SDSS 

Our goal is to identify pure samples of TV = 1 and TV = 2 
galaxies in order to investigate the differences in their prop- 
erties. We apply the selection criteria discussed in the pre- 
vious section to a volume-limited samples of galaxies from 
the New York University Value-Ad ded Galaxy Catalogue 
(NYU-VAGC, iBlanton et all |2005al ) based on SDSS DR6 
with M r — 51og 10 h < —19 (—20). The volume-limited sam- 
ples contain galaxies from the Main galaxy sample with an 
extinction-corrected apparent magnitude r < 17.64 in re- 
gions of redshift completeness greater than 0.8. The volume- 
limited sample with M r — 51og 10 h, <— 19 (— 20) contains 
67,472 (107,327) galaxies. 

First, we define a population of TV = 1 galaxies from 
the volume-limited sample. For each galaxy in the volume- 
limited sample, we identify all other galaxies within AV < 
1000 km s" 1 and compute Dn and TV7oo- We also compute 
the number of potential neighbours — galaxies in NYU- 
VAGC without a measured redshift but in the relevant mag- 
nitude range — within a projected distance of 700 h~ kpc. 
A population of TV = 1 galaxies is defined by selecting all 
galaxies with TV700 < 1 and Dn > 400/i _1 kpc and no po- 
tential neighbours. The SDSS -19 (-20) TV = 1 sample con- 
tains 7,266(26,192) galaxies. The expected purity of this 
population is /jv=i = 0.999 (0.995) with negligible amounts 
of contamination by TV = 2 and TV > 2 systems according to 
our mock volume- limited redshift survey (see Table [l}. 

A population of TV = 2 galaxies is defined by selecting 
all galaxies with TV700 = 1 and Dn < 200 ft -1 kpc and no 
potential neighbours. The SDSS -19 (-20) TV = 2 sample 
contains 918 (2,549) galaxies. The expected purity of this 
sample is /jv=2 = 0.709 (0.848) and the contamination is 
mostly from TV = 1 systems, /jv=i = 0.289 (0.144), according 



to our cosmological model. The contamination by systems 
with TV > 2, f N >2 = 0.002 (0.007), is very small and we 
will only directly consider and correct the contamination by 
TV = 1 systems. If instead we take the maximum Dn for the 
TV = 2 sample to be 250 or 150/i _1 kpc, the sample then 
contains 1,096 (3,229) and 721 (1,830) galaxies, respectively. 
However, the contamination is higher for 250 /i -1 kpc with 
/jv=i = 0.369 (0.203) and lower for 150 h' 1 kpc with f N=1 = 
0.199 (0.106). 

The overall TV = 2 population may be further divided 
into satellite and central galaxies. The satellite galaxy re- 
sides within a smaller subhalo, which has fallen into the 
larger host halo of the central galaxy. For each TV = 2 
pair, the less (more) luminous galaxy in M r is selected as 
the satellite (central) galaxy. 

3.2 Correction of TV = 2 systems for 
contamination by TV = 1 systems 

In the previous section, we observed that TV = 2 popula- 
tions selected using our constraints on Dn and TV700 are 
primarily contaminated by TV = 1 systems. Using our mock 
volume-limited redshift survey, we quantified the levels of 
contamination and the results are listed in Table Q] Below 
we discuss our technique to statistically correct an TV = 2 
population for contamination by TV = 1 systems. 

Let Hn=2{x) be the observed distribution of galaxy 
property x for an N — 2 population selected using con- 
straints on Dn and TV700. The fraction of TV = 2 systems in 
this population is /jv=2 and the fraction of TV = 1 systems is 
/jv = i, as estimated from the mock volume-limited redshift 
survey. Thus, 

H N =2(x) = /jv=i * H N u ll(x) + /jv=2 * Hfj u l r 2 (x), (1) 

where HfJ^Kx) is the distribution for a pure N — i popu- 
lation. We ignore the very small contamination by TV > 2 
systems. TV = 1 populations selected by our constraints on 
Dn and TV700 are >99% pure and H N xx ^(x) is known. The 
distribution for a pure TV = 2 population may be approxi- 
mated by 

ff pm/^ ^ H N=2 (x) - f N =i * H N Zl(x) , . 

- fJ JV=2l a; J 7 T~7 ■ \ z ) 

JN=2 + JN>2 

The satellite/central TV = 2 populations are also con- 
taminated by a population of TV = 1 systems. The popula- 
tions of TV = 1 systems contaminating the satellite/central 
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Figure 3. Contamination-corrected stellar mass distributions of the —19 (left) and the —20 (right) overall (top), satellite (middle), 
and central (bottom) N = 2 populations (dashed). Also shown is the stellar mass distribution of the full N = 1 population (solid) 
without Monte Carlo mass resampling. Without Monte Carlo mass resampling, differences in stellar mass contribute to difference in 
colour between N = 1 and N = 2 populations. 



© 2012 RAS, MNRAS 000,ITlfT5l 



Groups of two galaxies in SDSS: implications of colours on star fornation quenching 7 



N — 2 populations are not the same as the population con- 
taminating the overall N = 2 population. To apply Equation 
to the satellite/central N = 2 populations, the appropri- 
ate distributions, Hf^.i' cen (x) and HfJ 1 : L\' B!Lt (x), must be de- 
termined. When a system identified using our constraints on 
Dn and -/V700 is two N = 1 systems instead of an N = 2 pair, 
the less (more) luminous N = 1 galaxy is placed in the satel- 
lite (central) galaxy population. As a result, the N = 1 pop- 
ulation contaminating the satellite (central) N = 2 popula- 
tion is less (more) luminous than the population contaminat- 
ing the overall N = 2 population. We construct the N = 1 
populations contaminating the satellite and central N — 2 
populations from the N = 1 population contaminating the 
overall N — 2 population by first randomly choosing two 
galaxies within the full N = 1 population. The less (more) 
luminous galaxy is placed in the satellite (central) N = 2 
contaminant population. This process is repeated 10,000 
times. Figure [2] shows the stellar mass (top) and colour (bot- 
tom) distributions of the -19 (left) and -20 (right) N = 1 
populations contaminating the overall (solid blue), satellite 
(dashed green), and central (dot-dashed red) N = 2 popu- 
lations. 



3.3 Stellar mass distribution and resampling 

Here, we examine the stellar mass distribution of our N = 2 
and N = 1 galaxy populations and discuss our procedure 
for resampling these populations to be of the same stel- 
lar mass. The stellar mass distributions of the —19 (left) 
and —20 (right) overall (top), satellite (middle), and cen- 
tral (bottom) N — 2 populations (dashed) are shown in 
Figure [3] These distributions are contamination-corrected 
using Equation ((2J) and the distributions shown in Figure [21 
Also shown is the stellar mass distribution of the full N = 1 
population (solid). 

To determine the difference in colour due only to a dif- 
ference in star formation activity between two galaxy pop- 
ulations, we must compare populations of the same stellar 
mass. From the stellar mass distributions, we see that satel- 
lite N — 2 galaxies have less stellar mass when compared 
to N = 1 galaxies, on average. On the other hand, central 
N — 2 galaxies have more stellar mass when compared to 
N — 1 galaxies, on average. Galaxies w ith higher stellar 
mass content tend to be redder in colour (|Kauffmann et all 
2003). Thus, differences in the colour distributions of our 
N = 2 and N = 1 populations are in part due to differ- 
ences in stellar mass as well as differences in star formation 
activity. We devise a Monte Carlo technique to resample 
the much larger N = 1 population, i.e. randomly select- 
ing without replacement, subpopulations whose stellar mass 
distributions match the contamination-corrected N = 2 dis- 
tributions shown in Figure [3] Using our Monte Carlo re- 
sampling, we can ensure that we are comparing two pop- 
ulations of the same stellar mass, which we quantify using 
the Kologorov-Smirnov (K-S) test using the function Kolo- 
mogorovSmirnovTest in Mathematica, and the differ- 
ences in the colour distribution will only reflect differences 
in star formation activity. The K-S test probabilities when 
comparing the stellar mass of the —19 and —20 overall, satel- 
lite and central N = 2 populations to the full and resampled 
N = 2 populations are listed in Table [2] for one realisation 
of the Monte Carlo resampling. The resampled N = 1 pop- 



Table 2. Kolmogorov-Smirnov test probabilities when comparing 
the stellar mass of various N = 2 populations with the full and 
resampled N = 1 populations 





Mag 












Limit: 


-19 


-19 


-20 


-20 




N = 1: 


Full 


Resamp. 


Full 


Resamp. 


N = 


2 Overall 


0.002 


0.984 


0.000 


0.969 


N = 


2 Satellites 


0.000 


0.993 


0.000 


0.889 


N = 


2 Central 


0.000 


0.920 


0.000 


0.537 



ulations are clearly more like the various N = 2 populations 
in stellar mass than the full N = 1 population. 



4 RESULTS 

4.1 N = 2 colour distribution 

The goal of this paper is to investigate the differences in the 
colour between N = 1 and N — 2 galaxies and determine 
their implications on star formation activity. We compare 
the g — r colour of N — 2 populations with N = 1 popu- 
lations resampled to have the same stellar mass to ensure 
that colour differences are only due to differences in star 
formation activity. NYU-VAGC magnitudes and colours are 
standard Petrosian magnitudes and colours |Petrosianlll976l ; 
IStrauss et all [2 002) that are galactic extinction-corrected 
l|Schlegel et al.l " 19981) and _ff-corrected b y the template- 
fitting method of iBlanton fc Roweisl (|2007l ). 

The contamination-corrected g — r distributions of the 
— 19 (left) and —20 (right) overall N — 2 populations 
(dashed) and the resampled N = 1 populations (solid) of 
the same stellar mass are shown in the top row of Figure U 
We separate galaxies into red (non-star-forming) and blue 
(star- forming) based on a cut in g — r and measure the excess 
of red galaxies between the two populations. The difference 
between the red fractions of the N = 2 and N — 1 pop- 
ulations as a function of the g — r value used to separate 
the red sequence and blue cloud is shown in the bottom row 
of Figure [4] A g — r cut of 0.68 yields the maximum dif- 
ference between the red fractions. This point corresponds 
to the "green valley" and we define a red galaxy to be a 
galaxy with g — r > 0.68 and a blue galaxy to be a galaxy 
with g — r < 0.68. This cut is shown as a vertical black line. 
The red excess is simply the difference between the two red 
fractions. 

The g — r distributions of the N = 1 populations shown 
in Figure 2] change as the subpopulations selected during 
our Monte Carlo resampling changes. Thus, the red excess 
is different for each Monte Carlo realisation. We generate 
100 independent realisations of our Monte Carlo resampling 
technique and compute the red excess for each realisation. 
Taking the mean and standard deviation over all realisa- 
tions, we find that the overall N = 2 population has a red 
excess of 0.05 ±0.01 and 0.06 ±0.01 for the -19 and -20 
samples, respectively, relative to N = 1 populations of the 
same stellar mass. Thus, N — 2 galaxies are redder than 
N = 1 galaxies of the same stellar mass. The presence of 
just one other galaxy within the same dark matter halo is 
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Figure 4. (Top) Contamination-corrected g — r distribution of the —19 (left) and —20 (right) overall N = 2 population (dashed) and 
a resampled N = 1 population (solid) of the same stellar mass. (Bottom) Difference between red fractions of the —19 (left) and —20 
(right) N = 2 and N = 1 populations as a function of the g — r value used to separate the red sequence and blue cloud. Black vertical 
line shows our chosen red/blue separator value at g — r = 0.68 and the difference at this value is taken to be the red excess. 




Figure 5. (Top) Contamination-corrected g — r distribution of the —19 (left) and —20 (right) satellite N = 2 population (dashed) and 
a resampled N = 1 population (solid) of the same stellar mass. (Bottom) Difference between red fractions of the —19 (left) and —20 
(right) satellite N = 2 and JV = 1 populations as a function of the g — r value used to separate the red sequence and blue cloud. Black 
vertical line shows our chosen red/blue separator value at g — r = 0.68 and the difference at this value is taken to be the red excess. 
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enough to suppress star formation and induce redder colours 
to a limited extent. 

TV = 1 and N = 2 systems are the simplest environ- 
ments and one of the major advantages of studying such 
simple systems is that ACDM makes robust predictions for 
their histories. Tracking the merger histories in TV-body sim- 
ulations reveals that the differences in past major merger 
(<3:1) hi stories of TV = 1 a nd TV = 2 subhaloes are < 1 
per cent |Stewart et alj|2008h . Hence, mergers alone cannot 
account for the observed red excess and it is likely that an 
environmental process that quenches star formation is re- 
sponsible for the observed red excess. 

4.2 Satellite TV = 2 red excess 

We investigate the origins of the TV = 2 red excess of the pre- 
vious section by examining separately the satellite and cen- 
tral TV = 2 galaxy populations. The contamination-corrected 
g — r distributions of the —19 (left) and —20 (right) satel- 
lite TV = 2 populations (dashed) and the resampled TV = 1 
populations (solid) of the same stellar mass are shown in 
the top row of Figure [S] The difference in the red fractions 
of the satellite TV = 2 population and an TV = 1 population 
of same stellar mass as a function red/blue separator value 
is shown in the bottom row of Figure [S] As before we se- 
lect g — r = 0.68 as our red/blue separator value (shown by 
the black vertical line) and the red excess is the difference 
between the red fractions of the two populations. 

As before for the overall TV = 2 population, we gener- 
ate 100 independent realisations of our Monte Carlo resam- 
pling technique and compute the red excess of the satellite 
TV = 2 population relative to TV = 1 populations of the same 
stellar mass for each realisation. Taking the mean and stan- 
dard deviation over all realisations, which accounts for the 
variation in the red fraction of the TV = 1 population, we 
find that satellite TV = 2 populations have a red excess of 
0.15 ± 0.01 and 0.14 ± 0.01 for the -19 and -20 samples, re- 
spectively, relative to TV = 1 populations of the same stellar 
mass. In other words, 15 per cent of satellite TV = 2 galaxies 
have transitioned from the blue to the red sequence due to 
star formation quenching after being accreted into another 
halo, assuming TV = 1 galaxies are the progenitor of satellite 
TV = 2 galaxies in a statistical sense. 

Recent studies of satellite star formation quenching 
in galaxy groups reveal that quen ching must occur over 
long t ime-scales of order 2-3 Gyr jlKang fc van den Bosch! 

\, |We 



2008 ; Ivan den Bosch et al.l 2008 ; IWeinmann et all l2009f 



McGee et al.l 120091 . 1201 ll ). These authors suggest strangu- 



lation as the physical process at work in these systems. The 
hot halo gas of the infalling satellite galaxy is stripped (by 
ram-pressure or tides) and star formation slows as the cold 
gas is consumed. However, semi-analytic models show that 
instantaneous and complete removal of hot halo gas by ram- 
pressure stripping leads to a passive red fraction of satel- 
lite galaxies that is much h i gher than the observed frac- 
tion (IWeinmann et al. 2006bl : lKang"fc van den B osch 2008: 
iKimm et alJl2009h . This has lead some authors to simply de- 
crea se the stripping efficiency in their semi -analytic model 
fe.g. lFont et al.l2008l ; IWeinmann et al.l2010h or suggest tidal 



stripping of hot halo gas l|Weinmann et al.l l2010) to better 
match the observations. 

However, strangulation of hot halo gas by any mech- 



anism tends to produces too many galaxies in the green 
valley l|Weinmann et al"1l2010l ). Because of this I Wetzel et all 
(2012) do not support strangulation and point out that it is 
not clear that strangulation is efficient in low-mass haloes, 
such as the groups of two being studied here, as such haloes 
are not expected to have virial shock fronts which sup- 
port hot, virialized gas within the halo l|Dekel fc Birnboiml 
|2006J). They also find a persistent specific star formation rate 
(sSFR) bimodality, i.e. the lack of galaxies in the green valley 
at all halo masses. Taken together, they argue that the satel- 
lite quenching mechanism must bring about a rapid transi- 
tion from the blue cloud to the red sequence, which any form 
of hot halo gas strangulation struggles with. Ram-pressure 
stripping of cold gas is the natural quenching mechanism for 
a rapid blue to red transition. Although most estimates indi- 
cate that ram-pressure stripping is very inefficient in the low- 
mass haloes being studied here. lNichols fc Bland-Hawthornl 
l|201lf) suggests that cold disc gas can be puffed up by in- 
ternal star formation making the removal of cold gas by 
ram-pressure stripping possible even with a rarefied exter- 
nal medium of a low-mass halo. 

Following the semi-analytic approach, we include a sim- 
ple treatment of the removal of the cold gas from the in- 
falling satellite galaxies (possibly by ram-pressure stripping) 
to our cosmological model. First, we consider an immediate- 
rapid quenching scenario where cold gas is removed in- 
stantaneously and completely immediately upon accretion. 
This is modelled u si ng th e popula tion synth e sis m odels of 
iBruzual fc Charlotl (2003) with a IChabrierl (2003) initial 
mass function with model galaxies with exponentially de- 
caying star formation rates with r = 1 — 10 Gyr. The in- 
falling galaxy's star formation rate is set to zero when the 
galaxy reaches the "initial" g — r value and we record the 
amount of time it takes to reach to the red sequence for each 
r. Figure HJ] shows the average time to reach the red sequence 
(g — r > 0.68) after star formation has stopped as a func- 
tion of g — r at accretion. The data points and error bars 
represent the mean and standard deviation over all values 

Of T. 

Our analysis indicates that a galaxy with g — r ~ 0.3 will 
take an average of ~ 1.2 Gyr after accretion to reach the red 
sequence in the immediate-rapid quenching scenario. Any 
galaxy with g — r > 0.3 will take less time. From the means 
and standard deviations in Figure [6] we generate 100 realisa- 
tions of the transition time, transition (<7 — r), for all g — r bins 
between 0.3 and 0.68 by drawing from random distributions 
with the appropriate mean and standard deviation and then 
interpolating. 

Next, we determine the transition fraction, 
/transition (<? — r), for all g — t bins between 0.3 and 
0.68. In the immediate-rapid quenching scenario this is 
simply the fraction that have been within their host halo 
for at least the transition time computed above. The 
distribution of time spent in host halo for satellite galaxies 
in our cosmological model, f(t), is shown in the top panel 
of Figure [7] The bottom panel shows the right cumulative 
distribution, F(t) = X*t'>i f(t'), or the fraction that have 
been within their host halo for at least t Gyr. The transition 
fraction is given by /transition (g - r) = F (transition (5 - r)) 
and is computed for each g — r bin between 0.3 and 0.68 for 
each realisation. The satellite galaxies in our cosmological 
model have TV700 = 1, D N < 200 h' 1 kpc, and TV = 2. They 
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Figure 6. Time required to reach the red sequence after star for- 
mation stops (instantaneous and complete re moval of cold gas) ac- 
cordin g to the population synthesis models of lBruzual fe Chariot! 
(2003) using model galaxies with exponentially decaying star for- 
mation rates with r =1-10 Gyr. Data points and error bars rep- 
resent the mean and standard deviation over all values of r. 



are distinguished from central galaxies by a host/satellite 
flag tabulated in the catalog. We compute the look-back 
time at accretion for each satellite halo from the scale factor 
of the universe at accretion using fi m = 0.3, Qa = 0.7, 
h = 0.7. 

To predict the red excess, we require the transition frac- 
tion estimated above and the g — r distribution of a popula- 
tion of N = 1 galaxies that become satellite N — 2 galaxies. 
There is no reason to assume that N — 1 galaxies of a par- 
ticular colour are more likely to become a satellite N = 2 
galaxy. The g — r distribution of this population should be 
identical to the N = 1 population distribution shown in 
Figure [S] HN=i(g — r). The red excess is then given by 



red excess = (/] -Hjv=i(<7 — r) * /< 



transition ^ 



(3) 



where the summation is over the range 0.3 < g — r < 0.68 
and the average is over all 100 independent realisations. Our 
immediate-rapid quenching model predicts a red excess of 
0.62 ±0.01 for the -19 sample and 0.54 ±0.01 for the -20 
sample. The errors are taken to be the standard deviation 
over all 100 realisations. The red excess predictions of our 
immediate-rapid quenching model exceeded the value ob- 
served in the NYU-VAGC DR6 samples by a factor of ~ 4. 



While IWetzel et all l|2012l ) argue that the quenching 
mechanism must be rapid, it does not need to occur imme- 
diately upon accretion, i.e. a delayed-then-rapid scenario. 
We can easily model such a scenario by including a delay 
time before truncating star formation in our model galax- 
ies. The time to reach the red sequence (see Figure [6]) be- 
comes transition transition + tdelay As the transition times 
get longer, the transition fractions decrease and the com- 
puted red excess decreases. The red excess predicted from 
the model as a function of star formation quenching delay 
time is shown in FigureJS] A delay time of ~ 3.3 and 3.7 Gyr 
predicts a red excess in agreement with the observations for 
the —19 and —20 samples, respectively. 

A few caveats should be mentioned regarding the re- 
sults above. Our observational result indicates that present- 
day satellite N = 1 galaxies have a higher red fraction than 
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Time Since Accretion (Gyr) 

Figure 7. Time spent within host halo for simulated satellite TV = 
2 galaxies with Vmax > 196kms — 1 (solid) and V ma x > 140 km 
s _1 (dashed) in our hybrid JV-body/semi-analytic cosmological 
model. 
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Figure 8. Expected red excess of satellite N = 2 galaxies as a 
function of the delay time in star formation quenching when the 
satellite is accreted into another halo for a —19 (solid) and —20 
(dashed) mock redshift su rvey based on a model usi ng the popu- 
lation synthesis models of lBruzual fc Chariot! J2003T) and our hy- 
brid 7V-body/semi-analytic cosmological simulations. Also shown 
are the observed satellite N = 2 red excess for the —19 (dot- 
dashed) and —20 (dotted) SDSS samples. Error bars represent 
the lcr variation due to d ifferent star formation histories in the 
iBruzual fc Charlod l|2003h models. 
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Figure 9. (Top) Contamination-corrected g — r distribution of the —19 (left) and —20 (right) central N = 2 population (dashed) and 
a resampled N = 1 population (solid) of the same stellar mass. (Bottom) Difference between red fractions of the —19 (left) and —20 
(right) central TV = 2 and TV = 1 populations as a function of the g — r value used to separate the red sequence and blue cloud. Black 
vertical line shows our chosen red/blue separator value at g — r = 0.68 and the difference at this value is taken to be the red excess. 



present-day TV = 1 galaxies of the same stellar mass. In 
our immediate-rapid quenching model and delay-then-rapid 
quenching model, we use a present-day TV = 1 population to 
predict the red excess of present-day satellite N — 2 galaxies 
assuming the former are representative of the progenitors of 
the latter. Ideally, we would have used an TV = 1 popula- 
tion at the average redshift of accretion of satellite TV = 2 
galaxies instead of the present-day population. From Fig- 
ure [7] satellite N — 2 galaxies have been within their host 
haloes for 2.3 and 2.6 Gyr on average for the —19 and —20 
samples, respectively. The average redshift of SDSS galaxies 
is z ~ 0.1. Thus, we should be using an TV = 1 population 
at z ^0.3 2 and 0.34 for the —19 and —20 samples, r espec - 
tively (see iMcGee et al.ll201ll ). Ivan den Bosch et all (|2008l ) 
also compared two present-day galaxy population while as- 
suming one population was representative of the progenitors 
of the other. Following their arguments we do not expect sig- 
nificant evolution in the red fraction of TV = 1 galaxies at 
low redshifts and any evolution is probably towards lower 
red fractions at higher redshifts. As such, our model esti- 
mates are actually lower limits. 

Another assumption made in our simple model is that 
the orbits of infalling satellites are identical for satellites 
that are red or blue at accretion. If satellites that are red 
at accretion are biased towards particular orbits, then the 
distribution of time spent within the host halo may differ 
from that shown in Figure [7] and our red excess estimates 
will be affected. It is difficult to tell whether the red ex- 
cess estimates would be larger or smaller without "a priori" 
information on which orbits red satellites prefer. 

Finally, in this preliminary study we have treated 



gas removal as instantaneous and complete and we have 
not included AGN feedback in contrast to more so- 
ph isticated semi-analytic g alaxy formation models such 
as iKang fc van den Boschl (|2008l ). iFont et all (|2008h and 
IWeinmann et all l|2010l )~ Nevertheless, our delay-then-rapid 
model produces a star formation truncation time after ac- 
cretio n that is si milar to the 3 Gyr reported bv lMcGee et al.l 
Hon]) based on IFont et ail <|2008h . In future work, we will 
studies of the properties of groups of two galaxies in con- 
siderably more depth and address the short-comings of our 
simple semi-analytic model in order to gain a more thorough 
understanding of the process(es) at work in the simplest 
galaxy groups. 



4.3 Central N — 2 blue excess 

In the previous section we explored our satellite N — 2 
populations; we turn our attention to our central N — 2 
populations here. The contamination-corrected g — r distri- 
butions of the -19 (left) and -20 (right) central TV = 2 
populations (dashed) and the resampled TV = 1 populations 
(solid) of the same stellar mass are shown in the top row of 
Figure [9] The difference in the red fractions of the central 
N — 2 population and an TV = 1 population of the same 
stellar mass as a function of the g — r value used to sepa- 
rate the red sequence and blue cloud is shown in the bottom 
row of Figure [9] As before we select g — r — 0.68 as our 
blue/red separator value (shown by the black vertical line) 
and the red excess is the difference in the red fractions of 
the two populations. For central TV = 2 populations we find 
a blue excess of 0.06 ± 0.02 and 0.02 ± 0.01 for the -19 and 



© 2012 RAS, MNRAS OOO.ITHTol 



12 Trinh et al. 





0.8 1.0 



Figure 10. (Top) g — r distribution of a —19 (left) and —20 (right) satellite TV = 2 population (dashed) with a red central galaxy and 
a TV = 1 population (solid) of the same stellar mass from one realisation of our Monte Carlo resampling technique. (Bottom) Difference 
between red fractions of the —19 (left) and —20 (right) satellite TV = 2 population with a red central and an TV = 1 population of the 
same stellar mass as a function of the g — r value used to separate the red sequence and blue cloud. Black vertical line shows our chosen 
red/blue separator value at g — r = 0.68 and the difference at this value is taken to be the red excess. 





0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



9- 9~ 

Figure 11. (Top) g — r distribution of a —19 (left) and —20 (right) satellite TV = 2 population (dashed) with a blue central galaxy and 
a TV = 1 population (solid) of the same stellar mass from one realisation of our Monte Carlo resampling technique. (Bottom) Difference 
between red fractions of the —19 (left) and —20 (right) satellite TV = 2 population with a blue central and an TV = 1 population of the 
same stellar mass as a function of the g — r value used to separate the red sequence and blue cloud. Black vertical line shows our chosen 
red/blue separator value at g — r = 0.68 and the difference at this value is taken to be the red excess. 
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— 20 samples, respectively, relative to N = 1 populations of 
the same stellar mass from mean and standard deviation of 
100 independent realisations of our Monte Carlo resampling 
technique. Thus, for our chosen red/blue separator, central 
galaxies in a group of two are bluer than isolated galaxies of 
the same stellar mass. 

The red excess in satellite galaxies is likely due to star 
formation suppression from gas loss. The gas that is lost by 
the satellite may be fed directly to the central, triggering star 
formation. This transfer could explain the observed blue ex- 
cess in centrals. To investigate this hypothesis further, we 
examine satellite galaxies that are paired with red (blue) 
central galaxies and compare them with the overall satellite 
population. We begin by considering the stellar mass distri- 
butions of the three satellite populations: all satellites, satel- 
lites with a red central, and satellites with a blue central. 
These distributions are not contamination-corrected using 
Equation Q because we do not know the appropriate dis- 
tributions for correction. The three satellite populations are 
resampled to be of the same stellar mass using our Monte 
Carlo resampling technique and we generate 100 indepen- 
dent realisations. 

For each realisation, we compute the g — r distribution 
and red fraction for all three satellite populations. Figures 
H0l and [TT1 show the results for satellites with a red and blue 
central, respectively, compared to the overall satellite pop- 
ulation for one realisation. The red excess of the satellite 
population with a red (blue) central galaxy relative to the 
overall satellite population is computed for each realisation. 
From the mean and standard deviation over all realisations, 
satellites with a red central have a red excess of of 0.05 ± 0.03 
and 0.02 ±0.02 for the —19 and —20 samples, respectively, 
relative to the overall satellite population of the same stel- 
lar mass. Similarly, the population of satellites with a blue 
central have a blue excess of 0.00 ± 0.03 and 0.03 ± 0.02 for 
the —19 and —20 samples, respectively, relative to the over- 
all satellite population of the same stellar mass. In other 
words, red satellites slightly tend to pair with red centrals 
and blue satellites slig htly tend to pair with blu e centrals and 
"galactic conformity" (|Weinmann et al"]|2006al ) is somewhat 
present even in groups of two. However, we conclude from 
this analysis that direct gas exchange between the satellite 
and central is not a likely cause of the central blue excess. 

Triggered star forma t ion from a close sat ellite pass (e.g. 
iMihos fc Hernquisdll996l ; iBarton et al.ll2007D could also ex- 
plain the observed blue excess in centrals. To test this hy- 
pothesis, we divide the central N = 2 population into a 
sample with satellite separations < 100 h~ l kpc and another 
with satellite separations > 100 h~ kpc. This satellite sep- 
aration value divides the central population into two nearly 
equal sized subsamples. The subsample with the closer satel- 
lites contains a red excess of 0.08 ± 0.05 and blue excess of 
0.03 ±0.03 for the —19 and —20 samples, respectively, rela- 
tive to the subsample with wider satellite separations. Here, 
the errors are simply Poisson errors. We do not find that the 
central population with closer satellites to be bluer than the 
population with wider satellite separations. Hence, triggered 
star formation from a close pass is not a likely cause of the 
central blue excess either. 

Lastly, we note that the central blue excess could also 
result because galaxies with satellit es are also more l ikely 
to be actively accreting cold gas fsee lKeres et al.ll2009l . and 



references therein), which is usually associated with larger- 
scale filamentary over-densities in galaxy position. However, 
we have no means of directly testing this scenario presently 
and we do not attempt to substantiate any further. 



5 CONCLUSION 

In this paper, we have used a ACDM cosmological model of 
structure formation to devise a method for selecting galax- 
ies in that are isolated in their dark matter halo (N = 1 
system) and galaxies sharing their dark matter halo with 
exactly one neighbour (iV = 2 systems) based on Dn, the 
comoving projected distance to the nearest neighbour within 
AV < lOOOkms -1 and A^oo, the total number of galaxies 
within a comoving projected distance of 700 h~ kpc within 
AV < 1000 km s -1 . Our cosmological model enables us to 
understand and correct for the contamination by galaxies in 
other environments allowing us to study the full, uncontam- 
inated distributions of star-forming and morphological pa- 
rameters instead of just the average trends. Using a Monte 
Carlo resampling technique, we construct populations of iso- 
lated galaxies and groups of two galaxies with the same stel- 
lar mass distribution. Comparing the g — r colours of these 
populations gives a pure measurement of the difference in 
star formation activity between isolated galaxies and galax- 
ies in a group of two, which is the simplest group environ- 
ment and the simplest change in galactic environment. We 
find the following: 

(i) If galaxies are separated into a red sequence (g — r > 
0.68) and blue cloud (g — r < 0.68), N = 2 systems have a 
red excess of 0.05 ± 0.01 and 0.06 ± 0.01 for the -19 and -20 
samples, respectively, relative to N = 1 systems of the same 
stellar mass. Thus, the presence of just one other galaxy in 
the same dark matter halo is enough to quench star forma- 
tion to a limited extent. 

(ii) Examining the less luminous member, satellite N — 2 
galaxies have a red excess of 0.15 ±0.01 and 0.14 ±0.01 for 
the —19 and —20 samples, respectively, relative to N = 1 
galaxies of the same stellar mass. 

(iii) An immediate-rapid star formation quenching sce- 
nario, where the cold gas of the infalling satellite galaxy is 
instantaneously and completely removed immediately upon 
accretion yields a red excess prediction of 0.62 ± 0.01 and 
0.54 ± 0.01 for the —19 and —20 samples, respectively. Thus, 
an immediate-rapid star formation quenching scenario is in- 
consistent with the observations. 

(iv) A delayed-the n-rapid star format ion quenching sce- 
nario, as suggested bv lWetzel et alj l|2012h . with a delay time 
of ~3.3 and 3.7 Gyr for the —19 and —20 samples, respec- 
tively, yields a red excess prediction in agreement with the 
observed red excess for satellite N — 2 galaxies relative to 
N — 1 galaxies of the same stellar mass. 

(v) Examining the more luminous member, central N — 2 
galaxies have a blue excess of 0.06 ± 0.02 and 0.02 ± 0.01 for 
the —19 and —20 samples, respectively, relative to N = 1 
galaxies of the same stellar mass. 

(vi) Satellite N = 2 galaxies with a red central have a 
red excess of 0.05 ±0.03 and 0.02 ±0.02 for the -19 and 
—20 samples, respectively, relative to the overall satellite 
population of the same stellar mass. Satellite N — 2 galax- 
ies with a blue central have a blue excess of 0.00 ±0.03 and 
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0.03 ±0.02 for the —19 and —20 samples, respectively, rel- 
ative to the overall satellite population of the same stel- 
lar mass. Thus, red satellites slightly tend to pair with red 
centrals and blue satellites slightly tend to pair with blue 
centrals demonstrating that galactic conformity somewhat 
present even in groups of two galaxies. However, the central 
blue excess cannot be explained by a simple direct exchange 
of gas between the satellite and central. 

(vii) Central N = 2 galaxies whose satellite separation is 
< 100 hT 1 kpc have a red excess of 0.08 ± 0.05 and blue ex- 
cess of 0.03 ± 0.03 for the —19 and —20 samples, respectively, 
relative to central N — 2 galaxies with satellite separations 
> 100/i. _1 kpc. Thus, triggered star formation from a close 
pass is unlikely to cause the observed central blue excess. 
The central blue excess may be due to cold flows. 

Our most significant result in this preliminary study is 
that present-day satellite galaxies in a group of two have a 
higher red fraction compared to present-day isolated galax- 
ies of the same stellar mass. Thus, star formation quenching 
of satellite galaxies by the yet undetermined dominant group 
process occurs even in the simplest group environment where 
there are no other neighbours and galaxy-galaxy interactions 
with the central galaxy are unlikely. Taken together with 
our other results, we have demonstrated that environmental 
processes strongly influence even the sparest groups of lu- 
minous galaxies. Further investigation of the star formation 
rate, AGN fraction, age, metallicity and concentrations of 
these simple systems and direct comparisons with predic- 
tions from ACDM and high-resolution hydrodynamical sim- 
ulations will help undercover the causes of the trends seen 
here and contribute significantly towards a comprehensive 
understanding of galaxy evolution. 
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